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Abstract. - Work can be extracted from a single heat bath if additional information is available. 
For the paradigmatic case of a Brownian particle in a harmonic potential, whose position has been 
measured with finite precision, we determine the optimal protocol for manipulating the center and 
stiffness of the potential in order to maximize this work in a finite-time process. The bound on 
this work imposed by a generalized second law inequality involving information can be reached 
only if both position and stiffness of the potential are controlled and the process is quasistatic. 
Estimates on the power delivered by such an "information machine" operating cyclically follow 
from our analytical results. 



Introduction. — The idea of extracting work from a 
single heat bath has been investigated for a long time as a 
putative challenge to the validity of the second law of ther- 
modynamics. Thought experiments such as Maxwell's de- 
mon or Szilard's engine, which supposedly violate the sec- 
ond law, were resolved by taking into account the thermo- 
dynamic processing of the acquired information through 
Landauer's principle, which quantifies the cost of infor- 
mation erasure [11 [5]. Whether it is the acquisition or 
the deletion of this information that compensates for the 
gained work still generates debate [21 H]. The fact, how- 
ever, that measurements can help to extract work from a 
single thermal bath has been widely investigated, mainly 
in theory for the classical jSHS] as well as for the quantum 
regime [9Ul4]. and validated by a recent experiment |l5j . 
This idea is based on the strong link between statistical 
physics and information theory [TMT5] . in particular how 
the information-theoretic concepts of relative entropy or 
mutual information can be related to the thermodynamic 
notion of work . The generalization of the Jarzyn- 

ski relation [53] formulated in [21] and the subsequent re- 
finements of the fluctuation theorems [351 [55] take infor- 
mation explicitly into account as a way to extract more 
work than the corresponding free energy difference from a 
non-equilibrium process. 

In this Letter, we address the question of how to extract 
within a finite-time the maximum amount of work exploit- 
ing the information gained from one measurement, using 
the framework of stochastic thermodynamics f27l|28]. As 
a model system, we consider a Brownian particle in a har- 



monic potential as in an optical trap. Several groups have 
used such colloidal particles in optical tweezers for exper- 
iments illustrating and testing the concepts of stochastic 
thermodynamics [29U34j . We determine the optimal pro- 
tocol [35H39] for two cases where we either control only the 
center of a trap or additionally its stiffness. Only the lat- 
ter case allows a complete conversion of all the information 
gained by a position measurement into work. Finally, we 
estimate the power of cyclic engines based on such optimal 
protocols. 

Dynamics. — We study the dynamics of a one dimen- 
sional Brownian particle in a single heat bath of temper- 
ature T, trapped in a harmonic potential 



V{x,t) =k{T){x-\{T)f/2 



(1) 



and driven into non-equilibrium by controlling the time- 
dependence of the position A(r) and the stiffness fc(T) of 
this potential, from the initial values Ai = and fcj to the 
end values At and fcf during the total time t. We call sys- 
tem the particle in the harmonic potential and medium 
the surrounding heat bath. We neglect inertial effects 
such that the dynamics is governed by the overdamped 
Langevin equation 



X = -^id^V{x,T) + C(t), 



(2) 



the dot representing the time derivative and /i the mobility 
of the particle. The thermal white noise ^(r) with 



(C(r)) 
(C(T)C(r')) 



0, 



(3) 
(4) 
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obeys the usual fluctuation-dissipation relation, the brack- 
ets (...) representing the average of a quantity. 

Throughout the paper, we set four quantities to unity: 
the temperature T, the Boltzmann constant k^, the mo- 
bility ^, and the initial potential stiffness /cj. This corre- 
sponds to a scaling of energies by a factor k^T, entropies 
by fee, times by and lengths by ^Jk^T jki. This 

rescaling does not reduce the generality of our results, 
given that these four quantities are independent units. 

Since we want to analyze mean values of the work, we 
will not use the Langevin equation ([2]) but rather the cor- 
responding Fokker-Planck equation, which with our units 
reads 

drp{x, t) 9^ [dxV{x, t) p{x, t) + d^p{x, r)] 

= -dxj{x,T), (5) 

where p{x, r) is the probability distribution function 
(PDF) for the position of the particle and j{x, r) the prob- 
ability current. We further make the ansatz of a Gaussian 
PDF 



where Astot is the total entropy change and Asgys the 



p{x,t) 



1 



: exp 



^A4(6(r),2/2(r)) 



(6) 



with average 6(r) and variance y'^{T) defining the com- 
pact notation Afx{b,y'^). By inserting eqs. ([T]) and dU 
into the Fokker-Planck equation ^ and identifying the 
coefficients of the resulting polynomial in x, we obtain the 
two equations of motion. 



1 



fc(r) 

A(r) = 6(r) 



k{ry 



(7) 
(8) 



which relate the evolution of the control parameters A and 
k to those of the two moments of the distribution b and 
y. Eqs. (O and © allow us to transform the problem 
of finding the optimal protocols A*(t) and k*{T) into the 
one of determining the optimal functions 6*(t) and y*{T), 
which turns out to be much simpler. 

Thermodynamics. — On a single trajectory a;(r), 
the first law holds as 28i 



w[x{t)] = IW[x[t)] + q[x{T)], 



(9) 



where w > represents work applied in the system, which 
will either lead to an increase AV^ = V{x{t), t) — V{x{Q),Q) 
of the energy of the particle or be transmitted to the sur- 
rounding medium as heat q. In the following, we will 
talk of extracted work for w < 0. The maximum ex- 
tracted work thus corresponds to the minimum applied 
work. Following 40_ , we relate the heat g[a;(T)] to the 
entropy change Asm of the medium 

q[x{T)] = As„,[a;(r)] = As^oMt)] - As,y,[a;(T)] (10) 



change in the entropy of the system. By putting eq. ([T 
into the first law ([9]), we obtain 

w[x{t)] - AV[x{t)] + Astot[x(r)] - As,y,[a.(T)], (11) 

which expresses the work w applied along a trajectory x{t) 
as a function of the change of three quantities. We will 
minimize the mean total work 



W = {w[x{t)]) ~ / w[x{t)]p[x{t)] d[x{T)] 



(12) 



where the integration is over all possible trajectories x{t) 
with the associated probability p[x{t)]. We will not di- 
rectly compute ([T^ but instead use the fact that we can 
easily average the three terms on the right-hand side of 
cq. (ITT|) for time-dependent Gaussian distributions 

The mean energy change AE = {V{x, t)) — {V{x, 0)) of 
the particle after the process is given by 

= I Km - \,f + y'm] - h'm + r (o)i (u) 



with the initial conditions Ai = and ki — 1. The mean 
entropy of the system is the usual time-dependent Gibbs 



entropy Ssysir) 



[xir)]) 



J p{x, t) lnp(x, r) da;. 



The mean total entropy change of the system thus is 

A^,y, = ln(zj(i)/y(0)). (14) 
Once again following [40], we find that 



S'tot = (Stot) — 



p{x,t) 



dx^b'^ + y^ 



(15) 



so that, by integrating this rate over the duration t of the 
process, we obtain the mean change of the total entropy 



A^tot = / {b'+r)dT. 



(16) 



The mean total work follows by averaging eq. (|lip using 
eqs. (dSl), (O, and ^ as 



W 



^[(6(t)-Af)^ + y^(t)]-i[6^(0) + y2(0)] 



(17) 



In order to optimize expression (jl7p . we have to identify 
the initial distribution N'x{b{Q),y'^{0)). If the initial state 
is the thermal one, we have 6(0) = and y(0) = 1. How- 
ever, if we first perform a measurement, the initial distri- 
bution will depend on its outcome and allow us to drive our 
system such that work can be extracted. We will show that 
the second law still holds in the generalized form [^12,20] 



W > AF ~L 



(18) 
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where 



F = - In 



cxp {—V{x)) dx 



2 \ k 



(19) 



is the free energy of the system and / is the information 
acquired through the measurement as defined in the next 
section. 

Measurement. — The state of the system before the 
time T = of the measurement is represented by the PDF 
Pi{x) = Mx{bi, yf)- We assume that the measurement is in- 
stantaneous and that its outcome a;,,! is distributed around 
the true position x as 



p{x^\x) = JVx^{x,yl^), 



(20) 



where j/m is the precision of the measurement. From eq. 
([20|) we can extract the PDF of the measurement outcome 

p{xm) = J Pi(x)p(a;,„|x) da; = A/'^„(&i, yf + yf„). (21) 

Through Bayes' theorem, p{x\xm)p{xm) — Piix)p{xin\x) , 
we obtain the conditional PDF for the true position x of 
the particle for a measured 



p{x\xm) =J^x 



9 2 



Vt + ym yt + 2/n 



Xrayf + hyl 



(22) 



The distribution p(x\xm) is the initial distribution 
A/'k(&(0), ?/^(0)) that we have to use in eq. ([T7]) if wc have 
performed a measurement. 

A quantity of interest in our problem is the KuUback- 
Leibler distance or relative entropy (41j 



^(a;m) = / p(x|a;„i)ln 



■ In 1 



P{x\Xm 

p{x) 

yf 



2(yr + yl) V yf + m 



dx 



(23) 



between the distribution p(x\x^) and p{x). For a fixed 
measurement outcome Xm, this distance quantifies the dis- 
tinguishability of the two distributions. If we average eq. 
P51) over all possible measurement results x^ distributed 
with (|2ip . we obtain the mutual information 



p{Xni)I{Xm) dx 

n 



2 



1 + 4 



(24) 



where the bar denotes the average over Xm- The infor- 
mation that appears in eq. ([T^ is given by eq. (^5)) for 
a fixed measurement outcome Xm and by eq. (|24p if we 
average over Xm, respectively. 

Instantaneous processes. — We first illustrate the 
generalized second law (|18|) and show that by measuring 
once the position of the particle, we can extract work by 
moving the trap even though there is no free energy dif- 
ference involved. We instantaneously move the potential 



p(xm|x)/ 




p(x|xm) 



V(x,0*)\ 


1 /T\ 




y w 





X-fXm 



Fig. 1: Optimal instantaneous protocol. Left: Initially, the po- 
tential is centered at Ai = (black). The particle is in thermal 
equilibrium with its position x distributed according to Pcq{x) 
(blue, dashed). The measurement outcome Xm is distributed 
with p{xm\x) (red) with precision ym around the real position 
X. Right: The optimal protocol consists in instantaneously 
displacing the potential to the position Af , which is the center 
of the distribution p{x\xni) (green) and yields the work w, in 
average equal to H29p . 



from Aj = to the final position Xf{x^) according to the 
measurement outcome while keeping its stiffness constant 
{kf = ki — 1). In particular we want to find the optimal 
position Af (xni) which minimizes the applied work, i.e., 
maximizes the extracted work. 

For a system initially in thermal equilibrium, we have 
6i = and yi = 1 in the PDFs dUl) and ([22]). The work 
applied to the system is just the energy difference between 
final and initial state 

{x - Af(Xm))^ - X^ 



w{x,x,^) 



(25) 



For a fixed measurement outcome Xm, the mean work is 
given by 



W{x^) 



p{x\Xai)'w{x, Xni) dx 
Af(Xin) 



1 + yL 

which becomes minimal for 



ii + yf.r 



1 



yn 



(26) 



(27) 



This optimal value corresponds to the center of the dis- 
tribution p(x|xm), as shown in fig. [TJ The optimal work 
W*{xni) is obtained by putting eq. (P7)) into eq. (^5]) and 
is always negative, but still obeys the inequality (fT8)) since 
we have 



W*{x^)+I{x^) 



-1 



1 



2(i + y?n) 2 



In 1 



y 



>0 = AF. 
(28) 

This inequality holds for fixed Xm and then also averaged 
over Xm- Equality would be reached for y^ ^ oo, i.e., for 
an extremely imprecise measurement, for which we do not 
gain any work. If we average W*{x^) over all measure- 
ment outcomes x^ distributed with (|2ip . we obtain the 
mean optimal total work 

1 



W* 



p{x„i)W* (xa,) dx 

n 



2(1 + yi) 



(29) 
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For an error- free measurement, we would have j/m = and 
consequently W* — —1/2, which means that we were able 
to extract the total mean energy of the particle by apply- 
ing an instantaneous shift of the potential following this 
ideal measurement. However, the mean amount of infor- 
mation (j24l) acquired through the measurement is always 
strictly greater than the mean extracted work (I^Hl) . We 
will be able to convert more information into work only if 
we allocate a finite-time to the process. 

Optimal finite-time protocols. — We now investi- 
gate optimal processes occurring in a finite-time t and for 
which the end values Af and fcf of the control parameters 
are assumed to be fixed. The setup is similar to the one in 
[35] with the main difference that we first perform a mea- 
surement before driving the potential, which allows us to 
extract work from the system. 

In the following, the system is in thermal equilibrium 
before the measurement so that we have bi — and y[ = 1 
and consequently 



5(0) 



1 



Vn 



and 2/2(0) 



Vn 



1 



Vn 



(30) 



for the initial distribution (1221) after a measurement. We 
determine the optimal functions 6*(t) and 2/*(t) and 
through eqs. (O and ^ the optimal protocols A*(t) and 
fc*(T) that minimize the work (IT71) in two different cases. 

Case 1: Constant stiffness. We first analyze the case 
in which the stiffness ^(t) = fci = 1 is constant. The sub- 
script 1 will denote results obtained under this assump- 
tion. In this case, eq. ([7]) is simplified and we have 



y 



vy(o) 



(31) 



so that all terms depending on y in eq. ([T7|) cancel. The 
total work then assumes the form 



d7 



(32) 



that consists of one boundary term and one integral term 
with 6(0) given by eq. ([50]) . By optimizing first the in- 
tegral term with the Euler-Lagrange equations, we find 
that 5(r) must be a linear function. We therefore opti- 
mize eq. ((32|) with respect to the final value b{t) and find 
the expression 

26(0) + Xit 



2 + t 



(33) 



from which we can determine the optimal protocol XI{t) 
through cq. ([8|) and the optimal work Wi{xm) by insert- 
ing it into eq. (15^ . If we add this optimal work to the 
information (j23|) . we obtain 

yHO)-l , (6(0)-Af)2 



W^ixn,) + I{x^)^-liiyiO) + 



2 + t 



(34) 



which is strictly positive due to < y(0) < 1. The second 
law (dH]) is therefore respected since AF = 0. In every 
case though, the extracted work is strictly smaller than 
the gained information (j23p . even in the quasistatic limit. 
We will see that we can extract work corresponding to 
the full information if we also control the stiffness of the 
potential. 

Using eq. (|2T|) to average over all possible measurement 
results Xm, we find the mean total optimal work 



p{xm)Wi (x,^) dx„ 



2 + t 2{2 + t) I + yl 



(35) 



consisting of two terms, where the first one corresponds 
to the one obtained without measurement in [35) . The 
second one is always negative and corresponds to the work 
that can be gained through one measurement at constant 
stiffness. For i = 0, we have Wf{t = 0) = A^/2: there 
is only one way to go from A; = to a fixed Af, namely 
an instantaneous jump, and we do not have time to take 
advantage of the acquired information - in this case, the 
measurement is useless. In the quasistatic limit, t — oo, 
the first term vanishes and we can extract the work 

W(^^oo) = -7T7r^T^, (36) 



2(1 



Vn 



which is the same expression as the one found for instan- 
taneous processes in eq. (j29p , as there is no dissipation in 
both cases. However, the difference here is that we do not 
need to adapt the end position of the potential Af to the 
measurement outcome as it was the case in order to derive 
eq. (|29p . We can rather extract work for any fixed value 
of Af provided we allocate enough time t > 2(1 + yj^JAf to 
the process. 

For the particular value Af = 0, i.e., imposing a final 
return to the initial potential, eq. p5p shows that work 
can be extracted in average for any total time t of the 
protocol. In this case, the optimal work Wi is equal to 

62(0) 



t 



2 + t 2 



(37) 



which means that we extract work from any initial state 
and any outcome of the measurement. 

Case 2: Control of both parameters. If we control not 
only the position A(r) but also the stiffness A:(r) of the 
potential, we have to use the full expression p7)) for the 
mean work. Like we already did to optimize p2p . we can 
use the Euler-Lagrange equations to first optimize the in- 
tegral terms and show that 6(t) and j/(t) must be linear. 
The total work then reads 



Wo 



62(0) , (6(i)-6(0))2 



Af)2 



y^^m-mi_ 



In 



y(o) 



(38) 
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Fig. 2: Optimal finite-time protocols for ki — 1 and t — 1. 
(a) Optimal stiffness ^^(r) for three different values of (b) 
Optimal position A2(r) of the potential for j/m = 0.4, Af = 0, 
and three different measurement outcomes 



Here, 6(0) and y(0) are given by eq. ([50)1 and the subscript 
2 denotes the control of both parameters. We set Af = A; = 
and fcf = fci = 1 so that the final potential is equal to 
the initial one and AF — 0. The total work (I55|) can be 
separated in two terms, one depending only on b{t) and the 
other on y(t). Some arithmetic yields the optimal values 
for the final distribution. 



bm = 



26(0) 

2 + t' 

2 + t 



(39) 
(40) 



that we can put into eq. (|38p to obtain the optimal work 
and into eqs. (O and ^ to get the optimal protocols, 
shown in fig. O The protocols show two jumps at r = 
and T = t, as previously found in |35| . The first jump 
can be quantified by putting the moments pop of the ini- 
tial PDF into eqs. ([7]) and This first jump adapts 
the shape of the potential to the measured position of the 
particle, depending on the precision of the measure- 
ment - the more precise the measurement is (y,„ 0), the 
sharper the potential after the measurement has to be. 
The particle is then dragged back closer to the position 
Af = while the potential also flattens towards fcf = 1. 
The end jumps vanish in the quasistatic limit t ^ oo. 

If we put the optimal values (15^]) and PO]) into eq. 




WT(") 



Fig. 3: Optimal mean work Wj 2 as function of the duration of 
the process t for = 0.3. If the stiffness k{r) of the potential 
is constant, the quasistatic work is given by eq. (|36|l . By 
controlling both A(r) and fc(r), we can extract the work (|42|) . 



and add the information given by eq. 



W^iXm)+Iix,n) 



([231) . we obtain 



2 + t 



t 



(41) 



+ 



2/2 - 1 



lny*(i) > = AF 



which is always positive since we have < ?/2(^) ^ 1- The 
generalized second law (fT8|) is therefore obeyed, as it was 
also the case for a constant stiffness. However, equality, 
W2* (2^111) + lixm) — 0, is now reached in the quasistatic 
limit since y^it — )■ 00) = 1. We can thus convert all ac- 
quired information into work in the quasistatic limit. By 
averaging W2 (a^m) over Xm, we obtain the mean total work 



W2 ■ Its limit for t cx) becomes 



Wiit 



00 



■In 1 



(42) 



and is exactly equal to the opposite of the mutual infor- 
mation (j24[) that we acquired through the initial measure- 
ment. The mean total work W2 is always greater than 
as shown in Fig. [H 
A second interpretation provides further insight into 
this process. In terms of extracted work, it is equivalent to 
a process without measurement from an initial potential 
of stiffness ki — l/y^(0) to a potential with ki = 1. In this 
case, the mutual information ([M)) is equal to the free en- 
ergy difference between the final and the initial state which 
could be extracted through a reversible (quasistatic) driv- 
ing of the control parameters. Furthermore, the optimal 
protocol is exactly the same, with the exception of the 
first jump, if we go without measurement from the initial 
potential 



Kn(x, 0) = - In ^2Try^0)p{x\x 



(x-6(0))^ 



(43) 



22/2(0) 

to the final potential Vf{x) — x'^/2. Thus, the measure- 
ment saves us some cost of preparing the system in the 
initial state (|43|) since this initial potential is adapted to 
the measured position Xm. It then becomes possible to 
extract work from the system for any duration t of the 
protocol. 
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Cyclic machines. — So far, we have analyzed the 
optimal protocol following one measurement. In a next 
step, it is interesting to conceive a periodically working 
engine which by repeated measurements extracts work. 
An analytical calculation of the power delivered by such 
a machine is slightly non-trivial because the PDF for the 
true (and hence the measured) position of the particle will 
depend on all previous measurements and protocols. A 
rough estimate, however, is possible using our results. If 
only the position of the potential can be controlled, eq. 
([36|) suggest that one can extract about the mean internal 
energy per relaxation time of the particle, i.e., a maxi- 
mum power P ~ ksTfiki if the units are restored. If in 
addition the stiffness can also be controlled, the maximum 
work per measurement according to eq. ()42p increases by 
a factor given by the logarithm of the precision of the mea- 
surement, which leads to the estimate for the maximum 
power P ~ kBTfiki\n{l/yU for < 1. 

Concluding perspective. — For the paradigmatic 
case of a colloidal particle bound in a harmonic potential, 
we have studied how to extract work from a single heat 
bath using the information gained from a position mea- 
surement with finite precision. This work obeys a gener- 
alized second law inequality relating extracted work, in- 
formation and free energies. Equality can only be reached 
if both center and stiffness of the potential are controlled 
in a quasistatic process. For finite-time processes the in- 
evitable losses can be minimized by using the calculated 
optimal protocols. 

Our results should be accessible experimentally by us- 
ing a colloidal particle in an optical trap, a setup which 
has previously been used to determine distributions for 
heat and work that enter the various fluctuation theorems 
in stochastic thermodynamics. On the theoretical side, it 
would be interesting to extend our study of the relation- 
ship between the concepts of stochastic thermodynamics 
and those of information theory in finite-time processes 
beyond the harmonic case to systems with anharmonic or 
interacting degrees of freedom. 
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